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Abstract. The one-parameter nonrelativistic supersymmetry of Mielnik [J. Math. Phys. 25, 
3387 (1984)] is applied to the simple supersymmetric model of Caticha [Phys. Rev. A 51, 
4264 (1995)] in the form used by Rosu [Phys. Rev. E 55, 2038 (1997)] for microtubules. By 
this means, we introduce Montroll double-well potentials with singularities that move along 
the positive or negative traveling direction depending on the sign of the free parameter of 
Mielnik's method. Possible interpretations of the singularity are either microtubule associated 
proteins (motors) or structural discontinuities in the arrangement of the tubulin molecules. 



PACS number(s): 87.15.He, 03.65.Ge, 11.30.Pb 

Microtubules (MTs) are hollow cylinder tubes, 25 nm in outer diameter and 14 nm inner 
diameter, made of two types of 4 nm-long dimers of a polar protein known as tubulin that 
can self-assemble both in vivo and in vitro to lengths from several nm up to mm in some 
neurons. They form the main filamentary component of the cytoskeleton of all eukaryotic 
cells. Along their walls the tubulin dimers are distributed onto 13 (the seventh Fibonacci 
number) so-called protonTaments laterally associated. Brain tissues are especially enriched in 
MTs. Many interesting speculations on MTs have been advanced in recent years Q- 

Based on well-established results of Collins, Blumen, Currie and Ross Q regarding 
the dynamics of domain walls in ferrodistortive materials, Tuszyhski and collaborators OJIU 
considered MTs to be ferrodistortive and studied kinks of the Montroll type [5 1 as excitations 
responsible for the energy transfer within this highly interesting biological context. 

The Euler-Lagrange dimensionless equation of motion of ferrodistortive domain walls as 
derived in [2| from a Ginzburg-Landau free energy with driven field and dissipation included 
is of the traveling reaction-diffusion type 



where the primes are derivatives with respect to a traveling coordinate £ = x — vt, p is a 
friction coefficient and a is related to the driven field 0. 

There may be ferrodistortive domain walls that can be identified with the Montroll kink 
solution of Eq. © 



ip" + pip' - ip 3 + ip + o = , 



(1) 




where (3 = (a 2 
cubic equation 



— ai)/v2 and the parameters ot\ and a 2 are two nonequal solutions of the 



(iP 



- a 2 )(ip - a 3 ) = ip 3 - ip - a . 



(3) 
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In a previous paper ©, one of the authors noted that Montroll's kink can be written as a 
typical tanh kink 



Af (£) = 7 - tanh 



(4) 



where 7 = 0:1 + 0:2 = 1 + 9a p- The latter relationship allows one to use a simple construction 
method of exactly soluble double-well potentials in the Schrodinger equation proposed by 
Caticha [7|. The scheme is a non-standard application of Witten's supersymmetric quantum 
mechanics JH1 having as the essential assumption the idea of considering the traveling M kink 
as the switching function between the two lowest eigenstates of the Schrodinger equation with 
a double-well potential. Thus 

0i = M0 O , (5) 

where O1 are solutions of (j>Q X + [e ,i — w(O]0o,i(O = 0, and tt(£) is the double-well 
potential to be found. Substituting Eq. © into the Schrodinger equation for the subscript 
1 and substracting the same equation multiplied by the switching function for the subscript 0, 
one obtains 



<f) + Rm4>o = , 
where Rm is given by 

M" + eM 



R 



M 



2M' 



(6) 



(7) 



and e = e 1 — e is the lowest energy splitting in the double-well Schrodinger equation. 
In addition, notice that Eq. © is the basic equation introducing the superpotential R in 
Witten's supersymmetric quantum mechanics, i.e., the Riccati solution. For Montroll's kink 
the corresponding Riccati solution reads 



Anh(K) 

and the ground-state Schrodinger function is found by means of Eq. © 



2 7 cosh 2 



(8) 



>0,M 



(0 = O (O) cosh (P-t) exp 



2f3 2 



exp 



2/3 2 



cosh(/3£)-7/3£-7sinh(/3£) 



(9) 



while is obtained by switching the ground-state wave function by means of M. This 
ground- state wave function is of supersymmetric type 



0o,m(O = 0o,A/(O)exp 



- / Rhi{y)dy 
Jo 



(10) 



where </>o,m(0) is a normalization constant. 

The Montroll double well potential is determined up to the additive constant e by the 
'bosonic' Riccati equation 



u M {i) = R M - Rm + eo = -r + 



d 2 , (7 2 -l)e 2 



+ 2 +eo+ 



4 ' 4/3 2 
(4 7 2 e + (2 7 2 + l)ecosh(/?0 - 8(3 2 ) cosh(^) 
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4 7 (e + ecosh(/3f ) - 2{3 2 ) sinh(/3£) 



(11) 



Plots of the asymmetric Montroll potential and ground state wavefunction are given in 
Figs. (1) and (2) for a particular set of the parameters. If, as suggested by Caticha, one 
chooses the ground state energy to be 



C ° 4 2 A(3 2 



(» 



7 



(12) 



then um(0 turns into a traveling, asymmetric Morse double- well potential of depths 
depending on the Montroll parameters (3 and 7 and the splitting e 



U, 



L,R 
0,m 



P 2 



1 ± 



2e7 



(13) 



where the subscript m stands for Morse and the superscripts L and R for left and right well, 
respectively. The difference in depth, the bias, is A m = Uq — Uq = 2e7, while the location 
of the potential minima on the traveling axis is at 



e, 



L,R 



1 . 



(2/3) 2 ±2e 7 
e(7Tl) 



(14) 



that shows that 7 7^ ±1. 

A one-parameter supersymmetric extension of the previous results is possible. It is quite 
known in the literature on supersymmetric quantum mechanics where it has been introduced 
by Mielnik, Fernandez and Nieto [9] and is based on the Darboux co variant isospectrality of 
Schroedinger equations. In the biological context it has been applied to the DNA molecule by 
Drigo-Filho and Ruggiero ifTUl . The point is that R M as given in Eq. © is only the particular 
solution of the Riccati equation ocurring in Eq. (ITTT) . A more general, parameter-dependent 
Riccati equation of the form A) = -Rf/(£; A) — R' M (£,; A) + e can be constructed whose 

solution is a one-parameter function of the form 



Rm(£; A) = R M (0 + ^[ln(/ A f (0 + A) 
and the corresponding one-parameter Montroll potential is given by 

u M fo A) = u M (0 - 2^ [ln(/ M (0 + A) 



(15) 



(16) 



In these formulas, Im{0 = J^0om(£)°^ an( ^ A is an integration constant that is used as 
a deforming parameter of the potential and is related to the irregular zero mode. The one- 
parameter Darboux-deformed ground state wavefunction can be shown to be 



0o,a/(£;A) = y / A(A+T) 



>0,M 



Im(0 + A 



(17) 



where ^A(A + 1) is the normalization factor implying that A ^ [0,-1]. Plots of A) and 
0o,m(£ ; A) for A = 10 are presented in Figs. (3) - (4). See also Fig. (5 ) for a plot of the function 
Im(0 producing the parametric Darboux deformation. A singularity at Im(0 + A = is 
introduced in both potential and wavefunction. If the parameter A is positive the singularity 
is to be found on the negative £ axis, while for negative A it is on the positive side. For 
large values of ±A the singularity moves towards =Foo and the potential and ground state 
wave function recover the shapes of the non-parametric potential and wavefunction. The one- 
parameter Morse case corresponds formally to the change of subscript M — > m in Eqs. (fl~5t 
and (fl"6l) . For the single well Morse potential the one-parameter procedure has been studied 
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by Filho [fTTl and Bentaiba et al [TT2"1. Potentials and wavefunctions with singularities are 
not so strange as it seems lfT3l . Similar to the case of the 5 potential in condensed matter 
physics, we interpret the singularity as representing the effect of an impurity moving along 
the microtubule in one direction or the other depending on the sign of the parameter A. In 
the case of microtubules, The impurity may represent a protein attached to the microtubule 
or a structural discontinuity in the arrangement of the tubulin molecules. This interpretation 
of impurities has been given by Trpisova and Tuszyhski in non-supersymmetric models of 
nonlinear microtubule excitations ifLH . 

In conclusion, the supersymmetric approaches allow for a number of interesting exact 
results and point to a direct connection between Schroedinger double-well potentials and 
nonlinear kinks encountered in nonequilibrium chemical processes. MTs are an important 
application but the procedures described here can be used in many other applications. 
Moreover, the supersymmetric constructions can be used as a background for clarifying 
further details of the exact models. Although it is not so clear why one should take a certain 
type of kink as switching function between the Schroedinger split modes, it is interesting that 
proceeding in this way one will be led to some familiar double-well potential in chemical 
physics. 
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Figure 3. The one parameter Darboux-modified MDWP for A = 10 



One-parameter nonrelativistic supersymmetry for microtubules 



6 




Figure 5. Plot of the integral Im(€) that produces the deformation of the potential and 
wavefunctions. 



